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Changes in the dynamics of supported polymer films in comparison to bulk materials involve a 
complex convolution of effects, such as substrate interactions, roughness and compliance, in addition 
to film thickness. We consider molecular dynamics simulations of substrate-supported, coarse¬ 
grained polymer films where these parameters are tuned separately to determine how each of these 
variables influence the molecular dynamics of thin polymer films. We Hnd that all these variables 
signihcantly influence the film dynamics, leading to a seemingly intractable degree of complexity 
in describing these changes. However, by considering how these constraining variables influence 
string-like collective motion within the film, we show that all our observations can be understood in 
a unified and quantitative way. More specifically, the string model for glass-forming liquids implies 
that the changes in the structural relaxation of these hlms are governed by the changes in the 
average length of string-like cooperative motions and this model is confirmed under all conditions 
considered in our simulations. Ultimately, these changes are parameterized in terms of just the 
activation enthalpy and entropy for molecular organization, which have predictable dependences on 
substrate properties and Him thickness, offering a promising approach for the rational design of film 
properties. 


I. INTRODUCTION 


Polymer films are used in a wide variety of applica¬ 
tions, ranging from micro-electronic devices to artificial 
tissues [1, 2]. However, both mechanical and dynamical 
properties of polymeric materials often change consid¬ 
erably in relation to bulk once confinement dimensions 
become less than ~ 100 nm. Much of the effort aimed 
at understanding the property changes in thin polymer 
films has centered on measurements related to the stiff¬ 
ness of these films [3, 4] and changes of molecular mo¬ 
bility, as quantified by the glass-transition temperature 
Tg [5, 6]. Many experimental [5, 7-14], as well as com¬ 
putational [15-20] studies, have reported large property 
changes in thin films. These changes have been mainly 
attributed to a combination of substrate interaction and 
geometrical confinement. There is also a growing aware¬ 
ness of the relevance of substrate roughness and stiff¬ 
ness, as well as non-equilibrium residual stress effects in 
cast films. It is a difficult matter to separate all of these 
different effects in experiments, and the present work ad¬ 
dresses this general problem through molecular dynamics 
simulations of substrate-supported, coarse-grained poly¬ 
mer melt films of variable thickness where the polymer- 
substrate interaction is varied, along with the boundary 
roughness and rigidity. Since we can tune these param¬ 
eters in simulations, we can obtain clear indications of 
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how each of these variables influence the film molecular 
dynamics. After an analysis of how these diverse factors 
affect basic dynamic properties of the polymer film, we 
show that the dynamical changes under all these condi¬ 
tions can be organized and understood in terms of how 
these constraining variables influence collective motion 
within the film, parameterized by the enthalpy and en¬ 
tropy of activation for molecular reorganization. 

Changes in Tg in polymer films are usually associated 
with local changes in the dynamics near the interfaces. 
Many studies have reported that a repulsive or neutral 
substrate along with a free boundary leads to an enhance¬ 
ment in dynamics and a reduction of Tg [8, 19, 21, 22]. In 
contrast, an attractive substrate, which typically slows 
down the dynamics near the substrate, results an in¬ 
crease in Tg [5, 8, 10-13, 15, 17, 20, 23-25]. However, 
an attractive smooth substrate with a relatively weak in¬ 
teraction may also enhance the rate of relaxation and 
diffusion [9, 11, 13, 15, 26, 27], demonstrating that the 
polymer-substrate substrate strength and the substrate 
roughness can also have significant effects on the poly¬ 
mer film dynamics. In particular, it has been noted that 
the enhancement or slowing of relaxation in supported 
films induced by two interfaces with different properties 
can complicate the interpretation of the thickness depen¬ 
dence of Tg [8, 11, 19-21, 23-25, 27]. 

The most prevalent type of polymer films are those 
supported on solid substrates, where the relaxation time 
is often increased near the substrate, while decreased at 
the free boundary. Additionally, experiments on multi¬ 
layered interfacial films have shown that the effects of the 
free boundary region can be largely eliminated by placing 
films between stacks of nano-layered polymers with differ¬ 
ent species [28] , suggesting that there is a length scale as¬ 
sociated with the interfacial film dynamics. This leads to 
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the question of whether the film dynamics depends sim¬ 
ply on the substrate interaction, or are there other phys¬ 
ically relevant characteristics of the interface that must 
be considered. After all, glass formation is a dynamical 
phenomenon, so that other variables - such as substrate 
rigidity - might be relevant. This motivates an explo¬ 
ration of the effects of substrate rigidity on properties of 
thin polymer films, a property that can be greatly tuned 
in polymeric materials through cross-linking or through 
control of the molecular structure [29-31] 

A popular picture to rationalize the changes of the film 
dynamics is a superposition of polymer layers with locally 
varying dynamics. In this simple model, any changes in 
the overall dynamics should be manifested locally. Thus, 
the interfacial layers are correspondingly expected to be 
the primary contributor to changes in the overall film 
dynamics. Near an attractive substrate, the polymers 
are ‘bound’ to the surface, leading to slower dynamics, 
while at the free boundary region of a supported or free¬ 
standing film, the chains have a relatively higher mobil¬ 
ity. At the film center, far from both interfaces, the local 
properties are expected to be ‘bulk-like’. This layer pic¬ 
ture of film dynamics is often conceptually linked to local 
changes in density profile or free volume. In our previous 
work, we found inconsistencies for this free-volume layer 
(FVL) rationale for explaining the observed changes in 
the dynamics [27]. Moreover, the dynamics can be en¬ 
hanced at the substrate, despite an increase in local den¬ 
sity. We also quantified the length scales of both density 
and dynamical perturbations within supported films and 
found that the temperature dependence of these scales 
are opposite to that at the free boundary region, incon¬ 
sistent with the FVL approach. The changes in the dy¬ 
namics of the film with a supporting layer are generally 
non-local, so it is naive to treat the film interior as being 
the same as bulk material. 

Here, we consider the dependence of the dynamics on 
film thickness, substrate roughness, and rigidity. We find 
that these parameters can induce significant changes in 
the dynamics, characterized by changes in the film Tg 
and fragility, but only rather subtle changes are observed 
in static properties, such as density. Again, we find that 
free volume ideas are not useful in predicting dynamics 
at the local level. Rather, substrate interaction, sub¬ 
strate roughness, and stiffness can all greatly influence 
the mobility gradient transverse to the substrate. Our 
findings for the variation of Tg with substrate rough¬ 
ness and interaction strength are consistent with earlier 
works [10, 15, 32, 33], but our observations on fragility 
and regarding substrate stiffness are new. Another novel 
aspect of the current work is that we characterize the 
fragility changes in terms of cooperative motion within 
the film and, in this way, obtain a quantitative under¬ 
standing of the wide variations in the temperature depen¬ 
dence of the structural relaxation time with film bound¬ 
ary conditions and thickness. 

There is continued interest in the breakdown of the 
Stoke-Einstein relation in glass-forming liquids and the 


possible relation of this phenomena with fragility and 
dynamical heterogeneity, and several recent studies have 
suggested specific relationships. Since we are able to 
tune the fragility over a large range using the same poly¬ 
mer model through modifications of confinement, we can 
asses the validity of these relations in our glass-forming 
polymer model. We find that the decoupling exponent 
relating the structural relaxation time to a diffusion re¬ 
laxation time can be systematically varied with confine¬ 
ment. The degree of decoupling increases as the effective 
dimension is reduced, i.e., smaller film thickness, consis¬ 
tent with recent observations for model glass-forming liq¬ 
uids in a variable spatial dimension [34]. Moreover, film 
fragility can either increase or decrease under confine¬ 
ment, depending on the boundary interaction, so we do 
not generally see an increase in decoupling with greater 
fragility, as suggested by crystallization measurements in 
non-polymeric materials [35]. Our results support recent 
observations [34] that indicate that changes in spatial 
dimensionality are relevant to understanding the decou¬ 
pling phenomenon. 

Given the sensitivity of the dynamics to the large col¬ 
lection of substrate properties, the question remains: how 
do we obtain a unified understanding of all these effects 
on the polymer dynamics? There has been much specu¬ 
lation that these changes revolve around changes in the 
collective dynamics of the polymer molecules, where the 
Adam-Gibbs theory is often discussed without a specific 
definition of the hypothetical ‘cooperatively rearranging 
regions’ (GRR) that are relevant to understanding these 
property changes. Simulations have identified cooper¬ 
ative rearrangements that are quantitatively linked to 
the structural relaxation time for bulk polymer mate¬ 
rials [36], and a similar connection has also been estab¬ 
lished in model polymer nanocomposites [37, 38]. These 
string-like motions therefore offer a molecular realiza¬ 
tion of the abstract GRR. We test this predictive scheme 
for the molecular dynamics simulations of polymer films 
where the inherent inhomogeneity of the dynamics of 
these materials makes it unclear whether is the model 
should still apply. Encouragingly, we obtain a remark¬ 
able reduction of all our simulation data for structural 
relaxation in thin polymer films based on this unifying 
framework. Lastly, we investigate the influence of con¬ 
finement on the activation free energy parameters that 
define our description. 


II. MODELING AND SIMULATION 

We model polymers as unentangled chains of beads 
linked by harmonic springs. The substrate is modeled ei¬ 
ther as a collection of substrate atoms, or by a perfectly 
smooth substrate. Non-bonded monomers or atoms of 
the substrate interact with each other via the Lennard- 
Jones (LJ) potential, and we use a shifted-force imple¬ 
mentation to ensure continuity of the potential and forces 
at the cutoff distance Tc- We choose Vc = 2.5aij to in- 
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elude inter-particle attractions where (Tij is the monomer 
“diameter” in the LJ potential. The index pair ij distin¬ 
guishes interactions between monomer-monomer (mm), 
substrate-monomer (sm), and substrate-substrate (ss) 
particles. The LJ interaction is not included for the 
nearest-neighbors along the chain. These monomers 
are connected by a harmonic spring potential C/bond = 
(r — ro) with bond length rg = 0.9 (equilibrium 
distance) and spring constant fcchain = (llll)emm/o'mm- 
rg. The spring constant is chosen as in Ref. [26], but 
we choose rg smaller than in Ref. [26] because we found 
crystallization occurs readily in the films for the value 
used in Ref. [26]. 

The interaction between monomers and the smooth 
substrate is given by. 


l^niooth 


27r 


riPs^s! 


15 




, ( 1 ) 


where 0 is the distance of a monomer from the substrate. 
This is the same smooth substrate model that we studied 
in our previous work [27]. To model the rough substrate, 
we tether the substrate atoms to the sites of triangular 
lattice (the 111 face of an FCC lattice) with harmonic 
potential. 



where rgq denotes an equilibrium position on the trian¬ 
gular lattice and kg is the harmonic spring constant [17]. 
We choose the lattice spacing to be 2^/®crss, where ags = 
O.SOCTnim and Usm = CTinm- All values are in reduced units, 
where a mm = 1 and emm = 1- Varying kg allows us to 
examine the role of substrate rigidity on the polymer dy¬ 
namics. We simulate films of variable thicknesses with 
Nc = 200, 300, 400, 600, or 1000 chains of 10 monomers 
each. These sizes correspond to thicknesses with value 
of roughly 6 to 25 monomer diameters. We use various 
interaction strengths (csm = s) between the rough sub¬ 
strate and polymers, ranging from 0.4 to 1.0 Cmm with a 
fixed surface rigidity kg = 100; we vary the strength of 
the substrate rigidity {kg = k) over the range from 10 to 
100 with a fixed e = 1. For this range of model param¬ 
eters, we find Tg of the film can be higher or lower than 
the bulk value. Additionally, we simulate a pure bulk 
system of 400 chains of M = 10 monomers each at zero 
pressure for the purpose of comparison. 

We define film thickness h{T) as a distance from the 
substrate where the density profile along the z direction, 
perpendicular to the substrate, p{z) decreases to 0.10. 
Other reasonable criteria does not affect our qualitative 
findings. The resulting h{T) is well described by an Ar¬ 
rhenius form, which we use to extrapolate the thickness 
value hg = h{Tg) at the glass transition. 

To quantify the overall dynamics of the films and bulk 
system, we evaluate the coherent intermediate scattering 


function. 


JV^( (3) 

where Vj is the position of monomer j and S{q) is the 
static structure factor. We define the characteristic time 
T by F{qQ,T) = 0.2, where qo is the location of the first 
peak in of S{q). To quantify dynamics locally within the 
film, we use the self (or incoherent) Fgeif (z, q, t) part {i.e. 
j = k) of Eq. (3) on the basis of the position z of a 
monomer at t = 0. We define the relaxation time t's(z) 
by Eseif(z,qg,rs) = 0.2. 


III. DEPENDENCE OF Tg AND FRAGILITY ON 
SUBSTRATE STRUCTURE 

A. Survey of Substrate Roughness and Film 
Thickness Effects 

We first contrast the overall changes to glass formation 
of polymer films with various thicknesses supported on a 
rough or smooth substrate having the same substrate- 
monomer interaction (e = 1.0). Relative to the bulk 
system, the relaxation time r of polymer films on the 
smooth substrate decreases as we decrease film thick¬ 
ness, and these deviations become more pronounced as 
we go to lower T, consistent with previous studies [26, 27] 
(Fig. 1(b)). However, we find the opposite trend for the 
rough substrate, as noted in Ref. [39, 40]. As we will see, 
this trend depends on substrate interaction strength and 
rigidity. We see that the dynamics change more rapidly 
with T for thinner films resulting in a larger r relative to 
the bulk material (Fig. 1(a)). We estimate Tg by fitting 
our data to the Vogel-Fulcher-Tammann (VFT) equa¬ 
tion. 


T(r) = Tooe^^«AT-To)^ (4) 


where Too is an empirical prefactor normally on the order 
of a molecular vibrational time (10“^'^ to 10“^^s) [41], 
D is a measure of ‘fragility’ and Tg is a temperature at 
which T extrapolates to infinity. Eq. 4 should only be 
applied above the glass transition temperature. In a lab 
setting, Tg is often defined as T at which the relaxation 
time reaches 100 s [42], and we adopt this simple criteria. 
Figure 1(c) shows that, relative to the bulk, Tg of polymer 
films on the rough substrate increases with decreasing 
film thickness, while for the smooth substrate systems, 
Tg decreases with decreasing film thickness. 

The variation in T dependence of relaxation is quan¬ 
tified by fragility, defined as the logarithmic slope of re¬ 
laxation time at Tg 


m{Tg) = 


cHnr 


d{TlTg) 


(5) 
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FIG. 1: Effects of film thickness and structure of the support¬ 
ing substrate on glass transition temperature Tg and fragility. 
The T dependence of relaxation time t of a bulk system 
and two representative film thicknesses hg supported on (a) 
a rough or (b) smooth substrate (e = 1). In this T range, it 
is apparent that, relative to the bulk, r increases as the film 
thickness is decreased on rough substrates, while for smooth 
substrates shows an opposite behavior, (c) Relative Tg and 
(d) relative fragility m to the bulk as function of film thick¬ 
ness. Both Tg and m for the rough substrates increase, while 
Tg and m of smooth substrates decrease as we decrease film 
thickness. 


We evaluate fragility m using the fit of Eq. (4). In 
Fig. 1(d), we see that, relative to the bulk, films on the 
rough substrate become more fragile as we decrease thick¬ 
ness, which is apparent from the increasingly rapid vari¬ 
ation of t{T) (Fig. 1(a)). In contrast, the fragility of 
polymer films on the smooth substrate decreases weakly 
with decreasing film thickness. 

Experimentally, Tg is often found to be proportional 
to m [43]. We also find a correlation between Tg and 
m for both substrates, but this relation is not strictly 
proportional. Note that films supported on a smooth 
substrate may have a non-monotonic thickness depen¬ 
dence of Tg and m on thickness. Specifically, our recent 
work [27] showed that Tg or m decreases with decreasing 
film thickness on the smooth substrate up to some crit¬ 
ical thickness, but that Tg increases for very thin films 
when interfacial effects become dominant. 


B. Local Structure and Dynamics 




FIG. 2: (a) Monomer density profile p(z) of a film, hg = 15, 
supported on a rough or smooth substrate, (b) Pair-pair 
correlation function in the direction parallel to the sub¬ 
strate (?(r||) near the substrate, (c) (/(r||) at the film center. 
Monomers near the rough substrate are slightly more densely 
packed and have better local ordering in comparison to those 
near the smooth substrate. 

To understand the observed changes in Tg and fragility, 
we resolve both structure and dynamics locally, since the 
changes in the properties of the film as a whole should be 
manifested in its local properties. We first contrast the 
local dynamics and monomer density as function of dis¬ 
tance z from the substrate boundary of rough or smooth 
substrates with monomer-substrate interaction strength 
e = 1 . We evaluate both p{z) and Ts{z) with a bin size 
6z = 0.875. 

In figure 2(a), we observe that the monomer den¬ 
sity near either the smooth or rough substrate increases 
weakly, and has a steady value through most of the film. 
The density drops to zero over a narrow window at the 
free boundary region. At the center of the film, the den¬ 
sity has a value close to the bulk. The density profile 
of the film on smooth substrate is essentially identical to 
that of a film on a rough substrate. 

In addition, we contrast the local structure parallel to 
the substrate by evaluating the density pair correlation 
function (/(r||) (see Fig. 2(b) and (c)). Far from the sub¬ 
strate, g{r\\) of both systems is indistinguishable, as the 
monomers are completely unperturbed by the substrate. 
Near the substrate, we see that there is a slight difference 
in the local structure. In particular, near the substrate, 
5 (r||) of the rough substrate has a somewhat larger first 
peak, indicating that the monomers near the rough sub- 
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strate are more ordered than those near the smooth sub¬ 
strate. In addition, there is a weak long-range ordering 
of monomers for the rough substrate, potentially induced 
by the periodicity of the substrate atoms. 



FIG. 3: Relaxation time Ts as function of distance a from 
the substrate Although the averaged densities of two systems 
are identical, the local dynamics is clearly distinct from one 
another, particularly near the the substrate. 

We next examine to what degree the local film dynam¬ 
ics reflect the changes in the density described above. 
Figure 3 shows that the dynamics of the film on a rough 
or smooth substrate at the same T are nearly identi¬ 
cal over the range from the center of the film to the 
free boundary region. However, there are large differ¬ 
ences of relaxation time near the substrate. The lo¬ 
cal relaxation time increases close to the rough sub¬ 
strate, but decreases near the smooth substrate. The 
enhanced dynamics near the smooth substrate are in 
part a consequence of the fact the monomers can “slide” 
along the substrate due to the substrate smoothness (see 
Refs. [27, 44]). This effect disappears for a rough sub¬ 
strate. An increasing relaxation time approaching the 
rough substrate has also been observed in a computa¬ 
tional study of a binary Lennard-Jones liquids, as well as 
in a bead-spring model of polymer melts with a relatively 
strong interaction, [16, 17, 39, 40]. Evidently, substrate 
roughness is highly relevant for the polymer film dynam¬ 
ics, and this factor must be controlled for consistent re¬ 
sults. 

A convenient way to parameterize local dynamical 
changes is by considering the local dependence of Tg and 
m as function of distance z from the substrate. This 
provides a way of summarizing the behavior of Ts{z,T), 
shown in Fig. 3. Figure 4 (a) shows that Tg increases 
near the rough substrate, reflecting the observed increase 
of Ts near the attractive substrate. Near the free bound¬ 
ary region, Tg decreases due to the enhanced mobility 
of monomers at the free boundary region. For relatively 
thick films, we find that there is a substantial film region 
where monomers have a Tg close to the bulk value. This 
is a situation in which the film thickness is large com¬ 
pared to the perturbing scales of the interfaces [27]. Tg 


is often found to be proportional to to, as observed in 
the overall dynamics. However, we do not see this pro¬ 
portionality between the local Tg and to. Specifically, m 
decreases approaching the rough substrate while Tg in¬ 
creases. This opposing trend has also been observed in 
polymer-nanoparticle composites [38]. 




z 


FIG. 4: Local Tg and fragility. (a) Relative Tg and (b) m 
of polymer films supported on rough substrate as function 
of distance z from the substrate for many thicknesses, (c) 
Relative Tg and (d) m of a polymer film {hg = 15) supported 
on a rough or a smooth substrate. Near the free substrate Tg 
decreases for both system. Near the substrate, Tg increases for 
film supported on rough substrate, but decreases for smooth 
substrate. 


Figures 4 (c) and (d) contrast the local variation of Tg 
and TO for rough and smooth substrates of a relatively 
thick film, hg = 15. In contrast to the increasing Tg 
of polymer films near the rough substrate, Tg of smooth 
substrate decreases close to the smooth substrate, which 
is consistent with variation of Tg (Fig. 3). Note that Tg 
and TO are slightly depressed for films on the smooth sub¬ 
strate, even at the middle of the film, a scenario where 
the perturbing scales of both interfaces become compa¬ 
rable to film thickness. 
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IV. DEPENDENCE OF DYNAMICS ON 
SUBSTRATE STRENGTH AND RIGIDITY 



2.5 1.5 2 

1/7 1/7 


FIG. 5: The T dependence of the relaxation time r for 
two him thicknesses, /ig = 15 and 6, with various interfacial 
strength e (a) and (b), and hxed rigidity k = 100. Panels (c) 
and (d) show the effect of variable rigidity k for hxed £ = 1. In 
general, r is decreased as we decrease the substrate strength 
or the molecular-substrate stiffness. 


A. Survey of Substrate Interaction and Rigidity 
Effects 

Substrate roughness is relevant to the film dynamics, 
but there are other crucial variables. We next investigate 
the dependence of dynamics on the interaction strength 
as well as rigidity of the rough substrate. First, we exam¬ 
ine the role of substrate interaction strength. Figures 5 
(a) and 5 (b) show how relaxation time r for two rep¬ 
resentative film thicknesses changes as we vary the inter¬ 
action e between the rough substrate and the polymers. 
The overall changes in dynamics result from the compet¬ 
ing effects of the substrate and free interface, so that r 
can be higher or lower relative to the bulk. As we have 
established, the free boundary region decreases r while a 
substrate with a relatively strong interaction increases r. 
Thus, for a given thickness, r decreases with decreasing 
the substrate interaction strength. 


We find a similar effect by varying the stiffness k of 
the bonds describing the substrate stiffness. Specifically, 
increasing the flexibility of the substrate (decreasing k) 
results in a smaller r (Fig. 5 (c) and (d)). Evidently, 
monomers of the chains near the substrate are less con¬ 
strained, since the substrate atoms are more flexible. The 
complete local analysis of the dependence of dynamics on 
flexibility of the substrate will be discussed in the next 
subsection. By comparing Figs. 5 (a) and (b), as well as 
(c) and (d), we can see that the substrate interaction or 
the flexibility of the substrate have greater influences on 
the thinner film, expected since the thinner film has a 
larger surface-to-volume ratio. 



E k 



E k 


FIG. 6: Dependence of the relative Tg and m of three repre¬ 
sentative film thicknesses as a function of substrate strength 
£ (a) and (b) or substrate rigidity k (c) and (d). For thin¬ 
ner films, the range of Tg and m is wide due to the larger 
substrate-to-volume ratio. 

We next evaluate the resulting dependence of Tg and 
fragility on the substrate interaction strength and rigid¬ 
ity of the rough substrate films. Figures 6 (a) and (b) 
show how Tg of three representative thicknesses change 
as a function of substrate interaction strength e at fixed 
rigidity k = 100. Generally, increasing the polymer- 
substrate interaction increases both Tg and fragility as 
monomer dynamics near the substrate presumably be¬ 
come progressively slower. These general trends of a de- 
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creasing Tg with decreasing substrate interaction have 
also been observed both in experiments and computa¬ 
tional works [10, 15, 32, 33]. This depression of fragility 
is also consistent with the findings in a free-standing 
film [18], which formally corresponds to taking the limit 
e —>■ 0. Evidently, the dependence of substrate polymer 
interaction of Tg or m becomes more significant for thin¬ 
ner films, as indicated by a steeper variation of Tg or m 
with s. 

We found similar trends for Tg and m by varying the 
substrate rigidity. That is, increasing substrate rigidity k 
at fixed substrate interaction strength (e = 1) increases 
both Tg and m. It is interesting to note that there ap¬ 
pears to be a nearly fixed point for Tg and m as func¬ 
tion of e. Specifically, Tg and m are independent of film 
thickness for e ~ 0.9 {k = 100) or fc ~ 75 (e = 1). We 
emphasize that this does not mean there are no changes 
in local dynamics, but rather that there is a balance be¬ 
tween the dynamic enhancement at the free boundary 
region and the slowing down of the dynamics near the 
substrate. In fact, the increasing behavior of m with de¬ 
creasing thickness is only observed for values k > 75 and 
e > 0.9. This compensation effect is reminiscent of the 
self-excluded volume interactions of polymers in solution 
near their 9 point [45, 46], and the compensation point 
for isolated polymers interacting with surfaces [47] . 

Both results potentially offer us insights into how Tg 
changes in multilayer films, which are “stacks” of poly¬ 
mer films with different species characterized by different 
flexibility, inter-polymer interaction, or molecular weight. 
Multilayer film experiments by Torkelson and co-workers 
have shown that a given layer of the multilayer film may 
have different Tg depending on the properties of neigh¬ 
boring layers [28]. Here, we emphasize that changes in 
dynamics do not necessarily arise from the substrate in¬ 
teraction strength alone; changes in the rigidity of the in¬ 
terface (e.g. polymer films placed on a polymer substrate 
with the same substrate interaction strength, but having 
different molecular flexibility) and substrate roughness 
are also relevant. 


B. Local Structure and Dynamics 

We revisit our analysis of both film structure and dy¬ 
namics (as in the previous subsection) to further confirm 
our arguments about the role of substrate changes on the 
overall dynamics. Figures 7 (a) and (c) show how the 
monomer density p{z) changes by varying the substrate 
strength or rigidity of the substrate. Far from either sub¬ 
strate, p{z) has a nearly constant value that is close to 
the bulk value. In general, the density near the substrate 
increases slightly, as we increase e or k. The behavior 
of the density at very low rigidity (k = 10) differs from 
more rigid substrates. These most flexible substrates can 
be thought of as an amorphous solid that does not per¬ 
turb the film much, because it can not adapt its structure 
to that of the polymer film. 






FIG. 7: Variation monomer of the density profile p(z) and lo¬ 
cal relaxation time Ts as function of distance z from the sub¬ 
strate with varying polymer-substrate interaction e (a) and 
(b) or substrate rigidity k (c) and (d). Note that p becomes 
nearly constant for 2 > 6, but for z < 6, p depends sensitively 
on the interaction at the boundary. 


Similar to our findings comparing rough and smooth 
substrates. Figures 7(b) and (d) show substantial changes 
in local relaxation Tg at the substrate as function of e or 
k. Local relaxation time Ts{z) generally decreases as we 
decrease e oi k. The weaker substrate interaction allows 
monomers to avoid caging near the attractive substrate. 
Likewise, decreasing substrate rigidity allows monomers 
to move freely, since the substrate atoms are not strongly 
localized. These dynamical changes do not mirror the 
changes in the local density. This again emphasizes the 
limitations of a free volume based interpretation of re¬ 
sults.. 


V. DECOUPLING AND THE ‘FRACTIONAL’ 
STOKES-EINSTEIN RELATION 

One of the canonical features of glass-forming liquids 
is that the decoupling of viscous and diffusive relaxation 
processes gives rise to a breakdown of the Stokes-Einstein 
(SE) relation approaching the glass transition. This ‘de¬ 
coupling’ phenomenon is frequently associated with the 






emergence of heterogeneity of the dynamics, which we 
know is prevalent in our thin polymer films. Normally, 
decoupling is quantihed by the relation between the dif¬ 
fusion coefficient D and viscosity or a collective relax¬ 
ation time. For polymer chains, D is not readily accessi¬ 
ble computationally since the mean-square displacements 
of the chain center of mass only reaches the diffu¬ 
sive regime after extremely long times when the polymer 
melt is cooled. Instead, reference [36] has offered evidence 
that the characteristic time t* at which the non-Gaussian 
parameter, 


a2{t) 


5(r2(t))2 


( 6 ) 


has a maximum provides a diffusive relaxation time that 
exhibits a decoupling relation to the a relaxation time r 
of the segmental dynamics. 

To broadly characterize the heterogeneity of segmental 
motion and estimate a diffusive time scale. Fig. 8 shows 
a 2 {t) for many T for hg = 15 with a rough or attractive 
substrate. Data for other thickness and substrate inter¬ 
actions show the same qualitative features, so we only 
show this representative example. As is widely appreci¬ 
ated, a 2 exhibits a peak at intermediate time t* . One un¬ 
usual feature of these data is that a 2 (t) does not decay to 
zero for large t, reflecting the fact that displacement per¬ 
pendicular to the substrate is intrinsically limited by film 
thickness. The inset of Fig. 8 shows the T dependence of 
t* for many film thicknesses for the rough substrate, and 
it appears that t* grows less rapidly on cooling for in¬ 
creasingly thin films. Again, data for different substrate 
interactions show the same general trends. 

We now use our data for t* and r to quantify the decou¬ 
pling of relaxation time scales. Typically, this decoupling 
gives rise to a ‘fractional Stokes-Einstein’ (fSE) relation 
described by a power scaling law, 

t* ~ (7) 


where (' < 1 is a fractional exponent characterizing the 
decoupling strength, so that C = 0 defines the simple case 
where the Stokes-Einstein relation is valid. Eigure 9 illus¬ 
trates this variation for the representative case of poly¬ 
mer films on a rough substrate with k = 100 and e = 1.0, 
where Q ranges from roughly 0.65 in the thinnest film to 
0.3 upon approaching the bulk limit. A similar increase 
of the decoupling strength C with increasing nanoparticle 
concentration has also been found in polymer nanocom¬ 
posites for both attractive and non-attractive polymer- 
nanoparticle interactions [48]. 

The decoupling phenomenon is an inherently more 
complicated problem in polymeric materials than for 
small molecule liquids because there are separate relax¬ 
ation timescales for the segmental motions within the 
chains and for center of mass motions (which occur at 
a much longer time scales associated with the displace¬ 
ment of the chain as whole). Recent work has shown 
that the fragility of the segmental and overall chain mo¬ 
tion relaxation processes are generally quite different [49] . 



FIG. 8: The non-Gaussian parameter 0:2 (t) for various T for 
the case of film thickness 15 with a rough, attractive substrate 
{e = 1, k = 100). The peak of 02 ( 1 ) defines the characteristic 
time scale t*. The inset shows the T dependence of t* for 
various film thicknesses. Data for other substrate interactions 
are qualitatively similar. 


Moreover, Ediger and coworkers [50] found a complete 
absence of decoupling between the center of mass diffu¬ 
sion and shear viscosity in unentangled polystyrene over 
a wide temperature range. This situation is contrasted 
with relaxation at a segmental timescale where we ob¬ 
serve a power law relation between t and t* . Sokolov and 
Schweitzer have separately considered a power law rela¬ 
tion between the segmental r and polymer chain relax¬ 
ation time, which they also described as being a ‘decou¬ 
pling’ relation [49]. We do not attempt to describe this 
result because the calculation of the chain relaxation time 
at low temperatures is computationally prohibitive, and 
because we do not believe that the ‘decoupling’ relation¬ 
ship of Sokolov and Schweitzer is analogous to the decou¬ 
pling relation found in small molecule liquids. Of course, 
this relation between the segmental and large scale chain 
dynamics is fascinating and deserves further study. 

We next consider experimental observations suggest¬ 
ing a direct relation between fragility and decoupling 
in small molecule glass formers. Decoupling in glass¬ 
forming liquids has mainly been studied in context of 
crystal growth [35, 51], where the decoupling exponent 
C is inferred indirectly from the relationship between the 
crystal growth cooling rate, and fluid viscosity, in par¬ 
ticular Ref. [35], shows a proportional relation between 
m and C for low molecular weight organic and inorganic 
glass-forming liquids for a wide variation in fragility [35] , 
although this work specihcally excluded polymeric mate¬ 
rials. Sokolov and Schweizer [49] studied the decoupling 
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FIG. 9: Parametric relation between the coherent relaxation 
time r and t* for various film thickness for the representative 
case of a rough substrate with k = 100 and e = 1.0, demon¬ 
strating a fractional power law relation t* ~ . Clearly, 

C, increases with decreasing film thickness as illustrated in 
Fig. 10. 


exponent relating the segmental relaxation time and col¬ 
lective chain motion relaxation time for different poly¬ 
mer glass forming materials based on dielectric relaxation 
time measurements, and found a monotonic increase C 
with fragility for this type of decoupling. Together, these 
works suggest that the degree of decoupling might gener¬ 
ally increase in systems having larger fragility, although 
there is no generally accepted theoretical understanding 
of why such a relation might exist. 

We compared this trend to our results in Fig. 10, where 
we observe that the decoupling strength increases with 
decreasing hg for smooth and rough substrates. The in¬ 
set of Fig. 10 illustrates how C varies with m, and we see 
that the relation between C and m does not follow a single 
trend in these two cases. A simple proportional relation¬ 
ship between C, and m does not describe our data, raising 
a question about the general relation between decoupling 
strength and fragility. However, decoupling seems to be 
uniformly enhanced by geometrical confinement. Sen- 
gupta et al. [34] have observed a diminished decoupling in 
small molecule liquids between D and r with an increase 
of spatial dimensionality. If we view making polymer 
films thinner as reducing the ‘effective’ spatial dimension, 
then our observations fully accord with those of Sengupta 
et al. Simulations of polymer nanocomposites also show 
a progressive reduction of fragility with nanoparticle con¬ 
centration [38] , an effect that might likewise be rational¬ 
ized by an effective dimensional reduction with increased 



FIG. 10: Dependence of the decoupling exponent on film 
thickness for rough and smooth substrate interactions when 
k = 100 and e = 1.0. For both cases, ^ increases on decreasing 
thickness regardless of the type snbstrate interaction. The 
inset shows the relationship between relative ^ and relative m 
values obtained by altering the thickness of the film. 

particle concentration, a point of view advocated previ¬ 
ously [52] . While a definite relation between fragility and 
the strength of decoupling seems unlikely from the find¬ 
ings of Fig. 10, dimensionality does seem to be relevant 
to understanding this effect. 

VI. COLLECTIVE MOTION IN THIN 
POLYMER FILMS 

It has long been argued that polymer relaxation is gov¬ 
erned by the scale of cooperative motion. From a theoret¬ 
ical perspective, the classical arguments of Adam-Gibbs 
(AG) [53] , and our extension of this model based on nu¬ 
merical simulation evidence and thermodynamic model¬ 
ing [54] , provides a theoretical perspective for testing this 
proposition. Specifically, according to AG theory, the ac¬ 
tivation Gibbs free energy AGa{T) is extensive of the size 
0 * of ‘cooperatively rearranging regions’ (GRR), so that r 
can be formally written in terms of the general transition 
state theory relation, 

r(T) = Toexp[AG,(T)/fcBr], AGJT) = z*Afi (8) 

where Afj, is the activation free energy at high temper¬ 
atures when particle motion does not involve a signif¬ 
icant cooperative motion so that z* equals a constant 
(AG originally assumed that z* ~ 1 at high temper¬ 
atures, corresponding to completely uncooperative mo¬ 
tion, but a constant value of z* at high T is all that is 
required to recover Arrhenius dynamics). Recent simula¬ 
tions have shown that, despite the rather heuristic nature 
of the original arguments of AG, Eq. 8 with z* identified 
specifically with the average size L of the cooperative 
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string-like particle exchange motion provides a good de¬ 
scription for t(T) in polymer melt simulations, even in 
the case when nanoparticles have been added to tune the 
fragility over a wide range [27, 36, 38, 55, 56]. Very re¬ 
cently, we have stablished a quantitative correspondence 
between the L(T) and a living polymerization theory [54], 
and inspired by these results. Freed [57] has systemati¬ 
cally derived Eq. 8 from transition state theory assuming 
that the transition states involve many-body transition 
events in the form of equilibrium polymers with z as the 
average string length. These results together provide a 
predictive theoretical framework for understanding the 
dynamics of glass-forming liquids. 

We next evaluate L{T) following methods described in 
previous works [36, 58, 59] to see if we can also describe 
the dynamics of thin polymer films within the same for¬ 
malism. Figures 11 (a) and (b) compare the L{T) for 
two film thicknesses with a rough or smooth substrate at 
two representative thicknesses. For a given film thick¬ 
ness, L(T) is larger and grows faster for a film with 
rough substrate, which is qualitatively consistent with 
r(T) (Fig. 1 (a) and (b)). The variation of L{T) with 
substrate interaction strength e or substrate rigidity fc, 
shown in Figs. 11 (c) and (d), and L{T) in these cases 
also qualitatively captures the variation of t{T) (Figs. 5 
(b) and (d)). The similarities between the variation of 
L{T) and t{T) suggest that we may be able to predict 
changes in fragility from the variation of L{T), as found 
in previous works [18, 37, 38, 55]. Accordingly, in the fol¬ 
lowing section we consider this possibility, and develop a 
framework to explain the variation of the activation pa¬ 
rameters relating L and r. 


VII. COLLECTIVE MOTIONS AS AN 
ORGANIZING PRINCIPLE FOR THIN FILMS 
DYNAMICS 


From the discussion of the preceding sections, it is 
clear that there are a variety of factors that can alter 
the dynamics of thin polymer films, including film thick¬ 
ness, roughness, the polymer-substrate interaction and 
the stiffness of the substrate. Impurities introduced from 
the film casting process and possible heterogeneity in the 
substrate chemistry due to, e.g., substrate oxidation, are 
also relevant [18]. These effects are all significant and 
the observed changes of the film dynamics involves the 
convolution of all these variables. We clearly need some 
organizing principle to explain how all these factors in¬ 
fluence the film dynamics, and guide the development of 
polymer films with rationally engineered properties. 

In this section, we explore a perspective that allows 
us to obtain a unified understanding of the diverse dy¬ 
namical changes. The Adam-Gibbs perspective of the 
dynamics of glassy materials generally emphasizes the 
importance of collective molecular motion to understand 
rates of structural relaxation, and the results of the pre¬ 
vious section indicate the promise of such an approach. 


In order to quantitively test Eq. 8, we follow Refs. [36- 
38, 60], and identify CRR size z* with the relative size 
LjLj^ of string-like cooperative particle arrangements. 
La = L{Ta) is the value of the string size at the tempera¬ 
ture Ta, above which an Arrhenius law for t(T) holds. To 
determine Ta, we use the same definition as in Ref. [54]. 

The analysis of the dynamics of our thin polymer films 
is then based on the string model for the dynamics of 
glass-forming liquids [54, 61], in which r is described by 
the AG inspired relation. 


r{h,T) 


To{h) exp 


L{T) Afi{h,T) 
La{^) ksT 


(9) 


where A^(/i,T) is the high temperature activation free 
energy for T > Ta, 


Afi{h, T) = AHa{h) - TASa{h), (10) 

where AHa{h) and ASa{h) are the enthalpic and the 
entropic contributions of the high T activation free en¬ 
ergy respectively. These basic energetic parameters vary 



FIG. 11: Panels (a) and (b) show the comparison of L{T) for 
films supported on a smooth or rough substrate. Panels (c) 
and (d) show the variation of tL{T) with varying polymer- 
substrate interaction e or substrate rigidity k. The variation 
of L{T) with hg, e and k mimics that of r(T) shown in Fig¬ 
ures 1 and 5. This qualitative consistency suggests the ap¬ 
plicability of the string model of relaxation [54] to quantify 
r{T). 
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with film thickness hg and type of interaction (e and 
k). Note that Eq. 9 for T > Ta_ becomes, T(h,T) = 
To exp[(A/i(ft,, T)/fcBT], the typical activation form of 
transition state theory. In fact, Eq. 9 at Ta implies that 
To(/i) is not a free parameter, but instead is determined 
by 


ro(/i) = T^(/i)exp[-A^(/i,TA(/i))/fcBr^(/i)], (11) 


where ta = t{Ta) so that AiJ^ and ASa are the only 
undetermined parameters in Eq. 9. This relation was 
noted and tested in Ref. [61]. The string model prediction 
for the structural relaxation time of a film of thickness h 
can then be formally written. 


r{h,T) 


rA(/i)exp 


L(T) A^,[h,T) 
LA{h) ksT 


Afi{h,TA) 

ksTA 


( 12 ) 

where, AHaih) and ASa{h) are the only parameters on 
which T depends, just as in ordinary transition state the¬ 
ory for homogeneous fluids. 



FIG. 12: Structural relaxation time r in terms of the aver¬ 
age strings size L for (a) various thicknesses, and (b) vari¬ 
ous polymer-substrate interactions or substrate rigidities, r 
is scaled by ro = ta exp[—A/r(rA)/fesrA] where A^{Ta) = 
AH a — TAASa, and AHa and ASa are determined by fitting 
to Eq. 9 over a broad T range. 


thus far studied. Figure 12 shows the linear relationship 
between ln(T) and A^L/ksT for different films thick¬ 
ness in panel (a) and for different substrate rigidity and 
strength of interactions in panel (b). The universal col¬ 
lapse of r in terms of string size was noted recently in 
a brief communications [55], but the variation of the re¬ 
laxation time prefactor ro(/i) was considered as a free 
parameter in that work. We find that this data reduc¬ 
tion holds for all film thicknesses supported on a rough 
or smooth substrates, and applies as well as to the bulk 
polymer material. The same reduced variable descrip¬ 
tion describes a representative film supported on a rough 
substrate for various substrate interaction or support¬ 
ing substrate rigidity (see Fig. 12 (b)). Although the 
data reduction is identical between these figures, we sep¬ 
arate them for clarity. This remarkable data reduction 
shows that we can quantitatively describe the film dy¬ 
namics of all these films based on the string model rela¬ 
tion (Eq. 9), despite a wide range of dynamical changes 
due to film thickness, polymer-substrate interaction, or 
substrate rigidity. For instance note that, Figs. 11(c) and 
(d) show that the average extent of cooperative motion 
L does not significantly change with e or k, but that 
the structural relaxation time does change considerably. 
Therefore, the changes in r must result from the varia¬ 
tions of AHa and ASa, which we next discuss. 


VIII. THEORY FOR THE CONFINEMENT 
EFFECTS ON ASa AND AHa 

In order to develop a theoretical model of relaxation in 
polymer thin films (or polymers with molecular additives 
and nanocomposites), we must understand what controls 
the basic activation parameters AHa and ASa in the bulk 
polymer reference system. In general, the activation pa¬ 
rameters at high T should depend on film thickness, since 
all these properties depend on the film thermodynamic 
properties. In this section, we specifically confront the is¬ 
sue of how the parameters AHa and ASa depend on film 
thickness, boundary geometry, and interaction strength. 
We begin by considering the variation of these energetic 
parameters in the high temperature limit, where cooper¬ 
ative motion does not complicate our discussion. 


A. Transition State Theory 

Classical transition state theory [62-64] implies that 
the diffusion coefficient, the structural relaxation time, 
and shear viscosity can all be described by an Arrhenius 
T dependence (at high T, where relaxation is not co¬ 
operative); i.e., the structural relaxation time r can be 
expressed by the Arrhenius expression 

T = To exp[AG'a/fcBT], (13) 


We now demonstrate the applicability of Eq. 9 in quan¬ 
titatively describing the dynamics of all the films we have 


where tq is the vibrational time, and the activation free 
energy AGa is associated with the displacement of a 
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polymer statistical segment [65, 66]. As discussed ear¬ 
lier, the activation free energy, AGa{T > Ta) = A/r, 
has enthalpic contributions AHa related to the strength 
of the intermolecular cohesive interactions, and entropic 
contributions ASa arising from entropy changes needed 
to surmount complex multidimensional potential energy 
barriers in condensed materials [67-70]. Predicting ASa 
is often a weak point in transition state modeling and 
the factor exp[—AS' q//cb] is often just absorbed into the 
measured prefactor tq as a practical matter, but this is 
not an option for glass-forming liquids. 

To guide our thinking, we need to recognize the phys¬ 
ical origin for the values of AH a and ASa- In order 
to understand qualitatively AH a, we go back to Eyring’s 
early transition state theory arguments [64] , and consider 
long-standing physical observations in simple fluids [71- 
73] that relate AHa to the energy change associated with 
the removal of a test molecule from its local environment 
in the fluid state [74-76]. Such, an interpretation has re¬ 
cently been implemented computationally by Egami and 
coworkers [77]. This perspective implies that AHa should 
scale in approximate proportion to the heat of vaporiza¬ 
tion i7vap or the cohesive interaction energy of the fluid. 
Although this argument is simple, its experimental va¬ 
lidity has been established for hundreds of fluids [78]. 
More recently, simulations of simple Lennard-Jones flu¬ 
ids in 2 and 3 dimensions have shown that AHa scales 
in proportion to the interaction parameter e [79-81], the 
natural measure of intermolecular interaction strength in 
simple pair potential models such as LJ fluids and also 
our polymer model. 

As noted above, the variation of ASa with molecu¬ 
lar parameters is less well understood. In many small 
molecule fluids, the intermolecular potential is weak, and 
therefore the variation of ASa can be reasonably ne¬ 
glected. However, for molecules with many internal de¬ 
grees of freedom, such as polymers, there can be a con¬ 
siderable variation in A^'^. In particular, a survey by 
Bondi [67] revealed that ASa/ks could vary over a 100 
units, and can even change sign, so that variations of ASa 
cannot be ignored. Bondi’s pioneering study describes 
the frustrations of early theoretical efforts to estimate 
ASa theoretically for complex fluids. 

The basic physical picture that the free energy of ac¬ 
tivation is related to the free energy cost of removing a 
molecule from its local molecular environment suggests 
a proportionate contribution to ASa from the cohesive 
intermolecular interaction [70, 82]. This effect is evident 
in Trouton’s rule [83, 84], which relates the heat of va¬ 
porization iJvap and the entropy of vaporization Svap of 
gases and the Barclay-Butler phenomenological relation 
linking enthalpies and entropies of solvation in many mix¬ 
tures [85-88]. Indeed, many studies have established the 
specific relation [68-70], 

A5„ = A^o + AHa/kBTaomp (14) 

supported by observations on diverse materials, where 
A^o captures a background contribution associated with 


the internal conhgurational degrees of freedom of fluid 
molecules. This linear relation has long been estab¬ 
lished for the Arrhenius activation parameters of bulk 
polymer fluids [89, 90] . In glass-forming materials, 
the entropy-enthalpy ‘compensation temperature’ Tcomp 
is often found to be near the glass transition tempera¬ 
ture of the fluid sometimes termed, the “melting tem¬ 
perature of the glass” [91]; in crystalline solid materials, 
Tcomp is often found to be near the melting tempera¬ 
ture Tjn of the solid [92, 93]. We indeed find ‘entropy- 
enthalpy’ compensation in our simulated glass-forming 
films (shown in Fig. 13) where the compensation tem¬ 
perature, Tcomp = 0.18, a temperature similar to the es¬ 
timated VET temperature. These observations suggest 
that Tcomp is determined by a physical condition at which 
the intermolecular cohesive interaction is insufficient to 
keep the material in the solid state, so that the fluid 
then begins to explore liquid-like configurations, but a 
quantitative understanding of how Tcomp relates to the 
structure of the potential energy substrate remains to 
be determined. We next examine how AHa and ASa of 
transition state theory depend on the scale of confine¬ 
ment in thin polymer films. 
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FIG. 13: Variation of AHa and ASa for different film thick¬ 
nesses, substrate roughnesses, substrate interactions, and sub¬ 
strate rigidities. The slope defines compensation temperature 
Tcomp = 0.18. The data for the thinnest film {hg = 6) shows 
a deviation from the thicker film data, which is discussed in 
the main text. 


B. Effect of confinement on AHa and ASa 

In general, film confinement (or the addition of ad¬ 
ditives to a fluid) must change the effective molecu¬ 
lar coordination number, and thus the cohesive interac¬ 
tion strength, so that ASa{h) and AHa{h) of activation 
also vary with film thickness and boundary interactions. 
This change in the effective coordination number due to 
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changes in the surface-to-volume ratio with confinement 
has been extensively discussed in connection to estimate 
critical point shifts in confined fluids and magnetic ma¬ 
terials [94, 95], and as the origin of the Gibbs-Thomson 
effect for the shift of the melting point of crystals [96]. 
This scaling argument suggests that AHa{h) scales with 
the substrate-to-volume ratio of the film, so that 

AHa{h) — Ai4a(bulk) ^ 1/h. (15) 

We see that this scaling is roughly consistent with the 
data in Fig. 14, but varying the substrate energy of the 
film can also lead to a change in AHa, whose sign depends 
on whether substrate-polymer interactions are stronger 
or weaker than the polymer-fluid interactions. In par¬ 
ticular, Fig. 14 shows that AHa decreases with k and e. 
Thus, there are a number of contributing substrate terms 
that influence AHa so a clean ^ 1/h finite size scaling 
of AHa is not obtained. Nonetheless, the general trend 
in variation of AHa with confinement is understandable. 
We emphasize that despite the highly variable nature of 
the variation of AHa and ASa with boundary interac¬ 
tion, roughness, substrate rigidity and film thickness, the 
compensation relation between these activation energies 
holds to an excellent approximation. However, there is a 
different relation between AHa and ASa for film thick¬ 
ness hg = 6, an effect that is also reflected of the increase 
in the decoupling exponent (/ in Fig. 9. This interesting 
effect requires further study, but it may be a consequence 
of the strong deviation of p from its bulk value very near 
the substrate. In particular, Fig. 7 shows large density 
gradients for z < 6, so that ‘ultrathin’ films with hg 
are rather unlike the bulk material and thicker polymer 
films. The component of ASa that is finked to AHa is 
expected from entropy-enthalpy compensation to follow 
the same dependence as AHa, suggesting a similar in¬ 
verse dependence on thickness. Fig. 15(c) shows this is 
the case. 

We have not previously discussed the relaxation time 
prefactor tq in Eq. 9, which varies appreciably with con¬ 
finement, as shown in Fig. 16 for thin polymer films on 
smooth or rough substrates. In particular. Fig. 16 indi¬ 
cates that the changes in tq can be as large as Id orders 
of magnitude, so this factor is highly relevant for un¬ 
derstanding the changes in relaxation time in thin films 
and nanocomposites. We again emphasize that tq is not 
a free parameter in the string model of glass-formation, 
but this quantity is entirely determined from AHa and 
ASa (Eq. 11). 

Although To varies strongly with confinement, the re¬ 
laxation time at Ta, ta varies only weakly, so that tq 
changes are due almost entirely to changes of AHa{h) 
and ASa{h). This phenomenon has not been appreci¬ 
ated before. Evidently, the significant changes in the 
relaxation of glass-forming films and nanocomposites de¬ 
rive in large part from the high temperature activation 
parameters, which are typically thought not to have a 
direct relation to glass formation. 






FIG. 14: Changes in the enthalpic contribution AHa of the 
activation free energy, (a) Effect of varying the roughness by 
varying feg. (b) Effect of varying the strength of the substrate 
interaction by varying e. (c) Changing confinement on smooth 
or rough films. 

IX. CONCLUSIONS 

We have systematically explored factors to that alter 
the dynamics of thin supported polymer films- film thick¬ 
ness, substrate roughness, polymer-substrate interaction 
strength, and the rigidity of the supporting substrate. 
All these factors were found to be highly relevant to the 
dynamics of our simulated polymer films and their cou¬ 
pling makes an understanding of changes in polymer film 
dynamics in thin films rather complicated. Simple free 
volume ideas are inadequate to explain the significant 
changes that we observe. Control of boundary roughness 
and polymer-substrate interaction is evidently necessary 
to make polymer films with reproducible properties and 
the prediction of film properties based on computation 
will require the specification of many factors related to 
the film boundary conditions and structure. 

Despite the wide variation of film dynamics with 
boundary conditions and film thickness, we find that we 
can obtain a remarkably general characterization of the 
changes in the film dynamics for all film conditions using 
the string model of structural relaxation. In particular, 
we quantitatively describe the change in dynamics of sup- 
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FIG. 15: Changes in the entropic contribution ASa of the 
activation free energy, (a) Effect of varying the roughness by 
varying the parameter fcg. (b) Effect of varying the strength 
of the substrate interaction by varying the parameter e. (c) 
Changing conhnement on smooth or rough films. 


ported polymer films based on how the collective motion 
is perturbed in the him. These changes are ultimately 
parameterized in terms of the high temperature activa¬ 
tion free energy, leading to the almost paradoxical finding 
that glass formation is controlled by the fluid properties 
in the high temperature limit. 


Given the computational success of the string model, 
we are now faced with the problem of determining the 
extent of string-like collective motion in real materials. 
Future work must address how the extent of collective 
motion can be effectively estimated from direct measure¬ 
ment. Our recent measurements [54] suggest a direct 
relation between the string length and the interfacial mo¬ 
bility scale near the polymer-air-boundary of supported 
films and offers a promising method for estimating the 
string length. Moreover, recent work [97] also suggests 
that noise measurements might be effective for estimat¬ 
ing average string length, and we plan to pursue this 
possibility in the future. 



FIG. 16: Changes in tq as a function of inverse of film thick¬ 
ness relative to the value for the bulk system, as a function of 
the inverse of the thickness hg. The circles represent smooth 
substrate and the squares represent those for a rough sub¬ 
strate. 


As a secondary consideration, we examined how the 
large changes in fragility in our simulations were related 
to changes between a diffusive relaxation time and the de¬ 
coupling relation in the segmental relaxation time. We 
find that decoupling and fragility can change in oppos¬ 
ing directions, calling into question a general relation be¬ 
tween decoupling and fragility. However, the changes in 
the decoupling exponent C are in accord with simulation 
results [34] and theoretical arguments [52] for the effect 
of reduced dimensionality on decoupling. Thus, our find¬ 
ings do not support previous observations indicating a 
proportional change between fragility and the strength 
of decoupling in small molecule liquids, and we conclude 
that the relation between fragility and decoupling re¬ 
quires further study. Our observations also imply that 
it would be worthwhile to define the concept of effective 
dimensionality more precisely and check this perspective 
with changes of the decoupling strength. 
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